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Abstract. Let G be a non-linear function of a Gaussian process \Xt\tei with long-range 
dependence. The resulting process {G(Xt)}tez is not Gaussian when G is not linear. We con- 
sider random wavelet coefficients associated with {G(X t )}tei. and the corresponding wavelet 
scalogram which is the average of squares of wavelet coefficients over locations. We obtain 
the asymptotic behavior of the scalogram as the number of observations and scales tends to 
infinity, ft is known that when G is a Hermite polynomial of any order, then the limit is 
either the Gaussian or the Rosenblatt distribution, that is, the limit can be represented by 
a multiple Wiener-Ito integral of order one or two. We show, however, that there are large 
classes of functions G which yield a higher order Hermite distribution, that is, the limit can 
be represented by a a multiple Wiener-lto integral of order greater than two. This happens 
for example if G is a linear combination of a Hermite polynomial of order f and a Hermite 
polynomial of order q > 3. The limit in this case can be Gaussian but it can also be a 
Hermite distribution of order q — f > 2. This depends not only on the relation between the 
number of observations and the scale size but also on whether q is larger or smaller than 
a new critical index q* . The convergence of the wavelet scalogram is therefore significantly 
more complex than the usual one. 
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1. Introduction 

Denote by X = {Xt}tez a centered stationary Gaussian process with unit variance and 
spectral density /(A), A £ (—tt,tt). Such a stochastic process is said to have short memory 
or short-range dependence if /(A) is bounded around A = and long memory or long-range 
dependence if /(A) — > oo as A — > 0. We will suppose that {X t }t& has long memory with 
memory parameter < d < 1/2, that is, 

/(A)~|A|- 2d r(A)asA^0 (1) 

where /*(A) is a bounded spectral density which is continuous and positive at the origin. This 
hypothesis is semi-parametric in nature because the function /* plays the role of a "nuisance 
function". It is convenient to set 

/(A) = |l-e- iA r 2d T(A), A€(-7r,7T]. (2) 

Since the process X is defined only if f(X)d\ < oo, we need to require d < \. 
Consider now a process {Y t }ti=z, such that 

(A K Y) t = G(X t ), t€Z, (3) 

for K > 0, where (AY) 4 = Y t — Y t -i, {X t }t£Z is Gaussian with spectral density / satisfying ([2]) 
and where G is a function such that E[G(Xt)] = and E[G(A() 2 ] < oo. While the process 
{Yt}t£Z is not necessarily stationary, its K—th difference A Yj is stationary and is the output 
of a non-linear filter G with Gaussian input. 
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We shall study the asymptotic behavior of the wavelet scalogram of jl ^ j^py,, that is, the 
avera g e of squares of its wavel et co efficients . As s hown in iFlandrin 1992 ]. Abry and Veitch 
[1998], Veitch and Abry 19991 ] and Bardet [2000] in a parametric context, the normalized 
limi t of scalogram can be used to estimate the long memory exponent d defined in ([1]). See 
also lAbrv et al,l 201 ll ] for recent empirical studies. 

In the semi-parametric conte xt, the case where the function G is linear was firstly con- 
sidered in iMoulines et al.l 2007] and the c ase where G is a Hermite polynomial of arbitrary 
order was studied in IClausel et all l2011bll. The case where G(Xt) is the so-called "Rosen- 
blatt process" was studied bv lBardet and Tudor! [2010l ] and is somewhat analogous to the one 
where G is the second Hermite polynomial. Our goal is to show that for more complicated 
functions G, one can obtain new types of limits. 

We have referred to Hermite polynomials a number of times. This is because they form a 
basis for the space of functions G and thus appear naturally in our setting. Since the function 
G satisfies E[G(X)} = and E[G(X) 2 } < oo for X ~ Af(0,l), G{X) can be expanded in 
Hermite polynomials, that is, 



G(X) = £ ■ 



jH q (X) . 



(4) 



One sometimes refer to 
sum flU) is in L 2 (Q), 

and 



as an expansion in Wiener chaos. The convergence of the infinite 



c q =E[G(X)H q (X)] , q>l, 



(5) 



are the Hermite polynomials. These Hermite polynomials satisfy Hq(x) = l,Hi(x) = x,H%(x) - 
x 2 — 1 and one has 

E[H q (X)H q ,(X)]= [ H q {x)H q ,{x)^=e- x2 l 2 dx = q\l {q=q , } . 
Jr vi7r 

Observe that the expansion @ starts at q = 1, since 

Co = E[G(X)Hq(X)] = E[G(X)] = , (6) 

by assumption. Denote by qo > 1 the Hermite rank of G, namely the index of the first 
non— zero coefficient in the expansion @. Formally, go is such that 

q = mm{q > 1, c q ^ 0} . (7) 

One has then 



+oo c 2 

£ = E[G(X) 2 ] < oo . 
Q=go 



(8) 



We will focus on the wavelet coefficients of the sequence {Y t } t& % in (J3j) - Since {Y t } te % is 
random so will be its wavelet coefficients which we denote by {Wj^, j > 0, k € Z}, where j 
indicates the scale and k the location. These wavelet coefficients are defined by 

W^ k = Y,H^-t)Y t , (9) 

where 7j f oo as j f oo is a sequence of non-negative decimation factors applied at scale j, for 
example 7j = 2 J and hj is a filter whose properties are listed in Appendix [Bl We follow the 
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engineering convention where large values of j correspond to large scales. Our goal is to find 
the distribution of the empirical quadratic mean of these wavelet coefficients at large scales 
j — > oo, that is, the asymptotic behavior of the wavelet scalogram 

Tlj—l 

S ^ = i E w h , (io) 

3 k=0 

adequately centered and normalized as the scale j and the number of wavelets coefficients rij 
available at scale j both ten d to inf i nity. 

The reduction theorem of iTaqquI 19751 ] states that if G(Xt) is long-range dependent then 

the limit in the sense of finite-dimensional distributions of X/fc=i G(X k ) adequately normal- 
ized, depends on the first term in the Hermite expansion of G. In other words, there exist 
normalization factors a n — > oo as n — > oo such that 

[nt] [nt] 

- Y, G{X k ) and - £ C ^H qo (X k ) , 

have the same non-degenerate limit as n — > oo. 

We are interested here, however, in the asymptotic behavior of the wavelet scalog rani S nij j 
in (fT0|) . We want to find exponents a > and v > such that as rij and j tend to oo, 

{nf 7 r»S n . +ud+u , U eZ}, (11) 

tends, after centering, to a limit in the sense of the finite-dimensional distributions in the 
scale u. This is a necessary and important step in developing methods for estimating the 
underlying long memory parameter. 

The limit of the sequence S nj j will be related to the so-called Hermite process. The 
Hermite process is a self-similar stochastic process, with stationary increments and long range 
dependence. The Hermite process of order q lives in the qth Wiener chaos, that is, it can be 
written as an iterated multiple integral of order q with respect to white noise. We refer to 
Definition 12.11 below for the precise representation. 

We will see that, in the scalogram setting, the reduction theorem mentioned above does 
not hold anymore. For example if G(Xt) = H\(Xt) + H qi (Xt), qi > 3 then the Hermite rank 
is qo = 1. But the limit of the normalized scalogram is not necessarily the same as that of 
Hi(X t ) = X t . This is essentially due to the fact that the scalogram involves squares and, in 
addit ion, depends o n two parameters j and rij which both tend to oo. 

In lClausel etaD |2011bl ]. the case 

G(X t ) = H q (X t ), q>2, 

was studied and it was shown that in this case the limit is a Rosenblatt process (see Defi- 
nition [2J]) . In the present paper we study other classes of functions G for which different 
Hermite processes appear in the limit. For example, for the process 

G(X t ) = H x {X t ) + H qi (X t ), qi > 3 , 

considered above, the limit of (jlip may be either Gaussian, a Hermite process of order q\ — 1 
or a Rosenblatt process depending on the specific circumstances. We will show the existence 
of a critical index q\ and of critical exponents u, v' such that when q\ < q\ , then : 

• the limit is Gaussian if rij <C 7J, 

• the limit is a Hermite process of order q\ — 1 if 7^ <C rij -C 7? , 
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• the limit is a Rosenblatt process if jj <C rij , 

where dj <C bj means that a,j = o(bj) as j — > oo. 

We will also study interesting cases where the function G has a Hermite rank greater than 
two. 

The paper is organized as follows. Long range-dependence and the multidimensional 
wavelet scalogram are introduced in Section [2j The main results are stated and illustrated 
in Section O The chaos decomposition of the scalogram is given in Section [U The study of 
the leading terms is done in Sections [5] and [6l The proofs of the main theorems are given in 
Section [7] while Section [8] contains some technical lemmas. Basic facts about the Wiener chaos 
are gathered in Appendix [A] and Appendix [B] lists the assumptions on the wavelet filters. 



2. Long-range dependence and the multidimensional wavelet scalogram 

The Gaussian sequence X = {X t }tez with spectral density ([2]) is long-range dependent 
because d > and hence its spectrum explodes at A = 0. Whether { H q (X t )}tpz is also 
long-ra nge dependent depends on the respective values of q and d. We show in lClausel et al.l 



2011a ], that the spectral density of {H q {X t )}te_z behaves like |A| as A — > 0, where 

5 + (q) = max(<%) , 0) and 8(q) = qd - {q - l)/2 . (12) 

Hence d+(q) is the memory parameter of {H q (X t )} t £z- Therefore, since < d < 1/2, 
{H q (X t )}tez, Q > 1) is long-range dependent if and only if 

%)>0«(l/2)(l-l/<z)<d<l/2, (13) 

that is, d must be sufficiently close to 1/2. Specifically, for long-range dependence, 

q = l^d>0, q = 2^ d> 1/4, q = 3 d > 1/3, q = 4 d > 3/8 . (14) 

From another perspective, 

5(q) > 1 < q < 1/(1 -2d) , (15) 

and thus {H q (X t )}tez is short-range dependent if q > 1/(1 — 2d). 

We shall suppose that the Hermite rank of G is qo > 1, that is the expansion of G(Xt) 
starts at qo. We always assume that {H qo (X t )}t<=z has long memory, that is, 

q < 1/(1 - 2d) . (16) 

The condition (|16p , with qQ defined as t he Hermite rank ([TP , e nsures such that {A K Y} teZ = 
{G(X t )}t& is long-range dependent (see IClausel et al.l 201 laf ] . Lemma 4.1). We are mainly 
interested in the asymptotic behavior of the scalogram S nj j, defined by (fTUj) as rij — > oo 
(large sample behavior) and j — > oo (large scale behavior). More precisely, we will study the 
asymptotic behavior of the sequence 

S nj+u j +u = S nj+n j +u -E(S nj+uJ+u ) = -— {W?+u,k-®(W? +Utk )) , (17) 

k=0 



adequately normalized as j,rij — > oo. 

There are two perspectives. One can consider, as in lClausel et al. 2011a ] . that the wavelet 
coefficients Wj +Ui k are processes indexed by u taking a finite number of values. A second 
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perspective consists in replacing instead the filter hj in (|9|) by a multidimensional filter hgj,£ 
1, • • • , m and thus replacing Wj t k in © by 



w *J,k = E h e,j(^jk - t)Y t . 



We adopted this second perspective in Clausel et al.1 1201 lbll and we als o adop t it here since it 
allows us to compare our results to those obtained in lRoueff and Taqqul [20091 ] in the Gaussian 
case. 

We use bold faced symbols Wj & and hj to emphasize the multivariate setting and let 
h i = {^j> ^ = 1, • • • , m}, W jjk = {W e j jk , £ = 1, • • • ,m} , 

with 

w i,fc = E h j(^ k ~ M = E h ihi k ~ t)A- K G(X t ), j > 0, k € Z . (18) 
We then will study the asymptotic behavior of the sequence 



s ^ = ^.E( w h-n^l k }) , (19) 

J fc=0 

adequately normalized as j — > oo, where, by convention, in this paper, 

Wh = {Wl hk , I = !,■■■ ,m}. (20) 
The squared Euclidean norm of a. vector x — [^l? • • • 5 

] T will be denoted by |x| 2 = x\ + 

• • • + rr m . 

1 1 turns out that the asymptotic behavior of S n > j depends on how the subsequence of Her- 
mite coefficients c q , q > 1 which are non-vanishing is distributed. We denote this subsequence 
by {c qe }e & £ where £ is a sequence of consecutive integers starting at 0, 

£C {0,1,2,...}, (21) 

with same cardinality as the set of non- vanishing coefficients, and (qe)e^c is a (finite of infinite) 
increasing sequence of integers such that 

go = index of the first non-zero coefficient c q , 

qi = index of the {£ + l)th non-zero coefficient , £ > 1 . 



Examples 

1) If 



G(X t )=c 1 H 1 (X t ) + ^H 3 (X t ) , 



where c\ 7^ 0, C2 = 0, C3 7^ 0, c q = for q > 4, then go = 1, gi = 3 and £ = {0, 1}. 

2) If 

G(X t ) = %H 2 (X t ) + C ^H 3 (X t ) + %Ht(X t ) , 

where c\ = 0, C2 7^ 0, C3 7^ 0, C4 7^ 0, c q = for q > 5, then go = 2, gi = 3, q2 = 4 and 
£ = {0,1,2}. 

3) If 

00 

c„ 



g=l 

where c q 7^ for g > 1 then go = 1, gi = 2, . . . , and £ = {0, 1, 2, • • • }. 
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4) If 

G(X t ) = , 

<?o! 



G(X t ) = -^H qo (X t ) 



where c qo ^ and c q = for q ^ qo, then C = {0}. 

While Co is always equal to (see ([6])), the assumption ([7]) ensures that c qo ^ and hence 
that C always contains the index 0, so that C is never empty. In particular, we may write 

(A K Y) t = G(X t ) = J2^H qe (X t ), t€Z, (22) 

where, if C is infinite, the sum converges in the L 2 sense. 
We set 

/ = {£ G C : £ + 1 G £, %+1 - % = 1} , (23) 

J = {(*l,40 G £ 2 : < £ 2 , %1 / 1 and q h - q h > 2} , (24) 

that is, g£ and qe+i take consecutive values when £ E I and and qg 2 differ by two or more 
when (^1,^2) G J. The structure of these two sets is particulary important. The set / could 
be empty (there are no consecutive values of qt) or not empty. Then we set 

^ [ min(I) > , when / is not empty , ,^g\ 
I 00 :, when I is empty . 

When £q is finite (that is, I is not empty), qi is the smallest index q such that two Hermite 
coefficients c q , c g +i are non-zero. It will be involved in the normalization. We define, in 
addition, 

m = mm({£ G C, qi > 3}) > . (26) 

Thus 

Qmo i s the smallest index q such that c q is non— zero with q > 3. 
Examples 

1) If 

G(X t ) = c x H x {X t ) + ^H 2 (X t ) + ^H 4 (X t ) , 

where c x / 0, c 2 # 0, c 3 = 0, c 4 / 0, c g = for g > 5 then C = {0, 1,2}, / = {1}, 
£ = 1, mo = 4 and J = {(2, 4)}. 

2) If 

G(X t ) = §# 2 (*t) + |^3(^) + %H*(X t ) , 

where ci = 0, c 2 / 0, c 3 7^ 0, c 4 / 0, c q = for q > 5, then £ = {0, 1, 2}, I = {2, 3}, 
£ = 2, m = 3 and J = {(2, 4)}. 

3) If 

G{X t ) = c X Hx{X t ) = cyX t , 

where ci 7^ and c g = for q > 2, then £ = {0} and both / and J are empty. 

We are interested in the asymptotic behavior of the normalized scalogram S n j defined 
in (|19p . This behavior depends on the sets J and I. These sets affect both the rate of 
convergence and the limit distribution of the rescaled sequence. The limit (see Section will 
be expressed in terms of the Hermite processes which are defined as follows : 
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Definition 2.1. The Hermite process of order q and index 

(l/2)(l-l/g)<d<l/2, (27) 

is the continuous time process 

f" p i(«lH \-u q )t _ 1 

Z q , d {t)= / ■ ■ r-|«l • • • u q \~ d dW(ui) ■ ■ ■ dW(u q ), t E M. . (28) 

Jri H VUq) 

It is Gaussian and called Fractional Brownian Motion when q = 1 and < d < 1/2. is 
non Gaussian and called Rosenblatt process when q = 2 and 1/4 < <i < 1/2. The marginal 
distribution of Zq^(t) at t = 1 is called the Hermite distribution of index q. It is called a 
Rosenblatt distribution when q = 2. 

The multiple integral ([2"8"]1 is denned in Appendix [A] The symbol f^ q indicates that one 
does not integrate on the diagonal Ui = Uj, j ^ i. The integral is well-defined when (|27p 
holds or equivalently when, 

1 < q < 1/(1 - 2d) , 

because then it has finite I? norm. This process is self-similar with self-similarity parameter 

H = dq+l — q/2 = 5{q) + 1/2 € (1/2, 1), 

that is for all a > 0, \Z q/ i{a t)}tm and {a H Z qj d(t)}teR have the same finite dimensional 
distributions, see iTaaaul [l9T9f ] . 

3. Main results 

We shall now state the main results and discuss them. They are proved in the following 
sections. We start with the assumptions 

Assumptions A {"W^, j > 1, k G Z} are the multidimensional wavelet coefficients defined 
by ([18]) , where 

(i) {Xt}t£Z is a stationary Gaussian process with mean 0, variance 1 and spectral density 
/ satisfying ((2|). 

(ii) G is a real- valued function whose Hermite expansion ()4]) satisfies condition (I16p . 
namely qo < 1/(1 — 2d), and whose coefficients in Hermite expansion satisfy the 
following condition : for any A > 

Cq = 0(e~ Xq ) as q -> oo . (29) 



(iii) the wavelet filters (hj)j>i and their asymptotic Fourier transform satisfy (W-a) 



(W-c) with M vanishing moments. See details in Appendix [ 



We shall focus on the asymptotic behavior of the scalogram for two basic classes of functions 

G. 

• The first class involves functions G with Hermite rank greater or equal to 2 and with 
two consecutive terms in the Hermite expansion, both of which having long-range 
dependence. The result is stated in Theorem 13.11 

• The second class involves functions G with Hermite rank equal to 1 with no two 
consecutive terms with long-range dependence. The results are stated in Theorems l3.2l 
and 13.31 

Other classes are left for future work. 
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3.1. G has a Hermite rank greater or equal to 2. Consider functions G of the form 



C2 



G(x) = ^H 2 (x) + ... + -^H geo (x) + 



9*0+1 



(s) + 



where ci = 0. Some of the c q , q >2 may be zero as well. More precisely assume that 

90 > 2 , (30) 

that is, that the Hermite rank of G is 2 or more. Also assume that (a) there exists two 
consecutive terms and that (b) both are long range dependent. Assumption (a) implies that 
the set I in (|23|) is not empty. Since the index qg (see ([25]) ) of the first of these two consecutive 
terms could be qo > 2, we have qt > 2. The index of the second of these consecutive terms 
is qi + 1 > 3. Assumption (b) will be satisfied if this second term is long-range dependent, 
that is 

q £o + 1< 1/(1 - 2d) , (31) 
by (I15p . We note that this situation implies the following boundaries for the parameter d: 

1/3 < d < 1/2 , 



as indicated in (|2j 
Set 

v = 2q £o + 1 - 2q . (32) 
The following theorem provides the limit of (|19f) for two different cases, depending on 



whether the limit of n ■ 7J when j — > +oo is null or infinite. It involves K > defined 
in ([3]), go m ©i $(q) is defined in (fT2"T) . £o i n ([25]) . The integer M is the number of vanishing 
moments of the wavelet filters and appears in Appendix [Bj 

Theorem 3.1. Suppose that Assumptions A hold with M > K + 5(qo). Suppose moreover 
that the Hermite expansion of G satisfies $3Q\) and \31\) . 

Then two limits in distribution of the centered multidimensional scalogram S n j in ( fi9j) . 
suitably normalized, are possible. They involve the Hermite processes in Definition \2.1\ evalu- 
ated at time t = 1. The coefficients involve £o and the multidimensional deterministic vector 
L q; whose entries [L q (hg j00 )]e=i,- ,m are defined as 



L q (he 



\he )00 (ui H Vu q 



\ui H h u, 



\2K 



n 

i=i 



U; 



-2d 



dui • • • du. 



(33) 



which is finite for any q < 1/(1 — 2d). Then 
(a) Ifrij <C then, as j,nj — > oo, 



n i_ 2d7 -2( % o)+^s^ > . (_5 (27r)90 -i [/.(OJftL^-x] Z V (1) 



(?0 



(b ) If <C f/ien, as j, nj — 

(l-2d)/2 -(5(% )+<5( %0 +l)+2/r)- 

Remarks. 



oo, 



{ S 2 V Y 1 ' (270*0 



r(0)% +l/2 L ]Z M (1). 
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1. Using (|132p with M > K and a > 1/2, t he integral in ( 331) is finite for any positive 



integer q < 1/(1 
tions (1301) and 




Thus, under Condi- 



appearing in the limits of Cases (a) 



and |(b)| have finite entries. 



In case |(a)[ the limit is a deterministic vector times the non-Gaussian Rosenblatt 
random variable Z2 J d(X)-> that is, the Rosenblatt process Z^&if) defined in (f28l) and 
evaluated at time t = 1. In case [Qoj| the limit is a deterministic vector times the 
Gaussian random variable Zi^(l), that is, Fractional Brownian motion Zi^(t) defined 
in (j28[) and evaluated at time t = 1. 

In the case where rij ~ C07J as j — )■ 00 for some Co > 0, the scalogram is asymptoti- 
cally a linear combination of a Rosenblatt and a Gaussian variable. Indeed, using the 
results of Section El one can see that the scalogram is the sum of two terms having 
the same order, both converging in the L 2 sense respectively to a Rosenblatt and a 
Gaussian variable. 



Proof. This theorem is proved in Section 17.11 



□ 



In the framework of wavelet analysis as in lMoulines et al.l 20071 ] . we have jj = 2 J and the 



number n 



n 



of wavelet coefficients available at scale j, is related both to the number N 



of observations Y\,- ■ ■ ,Yj^ of the time series Y and to the length T of the support of the 
analyzing wavelet. More precisely, one has (see Moulines et al.l |2007l | for more details), 



n ; 



x 



,j - [2-3 (N - T + 1) - T + 1] = N + 0(1) , (34) 

where [x] denotes the integer part of x for any real x. Note that the assumption rij 
when j — > 00 is equivalent to N — > 00 faster than 2 3 . Moreover, for any v > 0, 

nj < 2 ju ^ 2- j N < 2 jv N < 2 j(ly+1) when N 

Examples. We now illustrate Theorem 13.11 through three examples : 



00 



(35) 



(i) G 

(ii) G 

(iii) G 



H qo with qo > 2. 



H qo + H qo+1 with g > 2, q + K 1/(1 - 2d). 

H qo + H qo+1 + H qi with q > 2, q + K 1/(1 - 2d) and with q x - (q + 1) > 2, 
that is, J = {(q + l,<?i)}. 

In all cases, the integer go denotes the Hermite rank of G. 

Let us elaborate on the conditions on d and the resulting limits for these examples. For 
simplicity, we assume that the scalogram S n . j is univariate. 



Example (i) When G = H qo with qo > 2, I and J are both empty. Since / is empty 
one can regard lo and consequently qg and v as infinity, which suggests that we are in 

oo. The asymptotic behavior of 



case 



(a) independently of the growths of rij versus jj as j 
the scalogram of this example is treated by Theorem 3.1 in IClausel et ail |2011bl ] under the 
condition qo < 1/(1 — 2d). Indeed, the obtained rate of convergence is the same as in case (a) 
of Theorem 13.11 and the limit is also Rosenblatt. This also corresponds to the limit obtained 
by Bardet and Tudo r in th e case where Y itself is the Rosenblatt process (see Theorem 4 
of lBardet and Tudorl \20ld!\ ). 



Example (ii) Suppose G = H qo + H qo+ i, with qo > 2 and go + 1 < 1/(1 — 2d). Then 
J is empty and / = {go}- The Hermite rank of G is go and thus coincides with qe . As a 
consequence, by (|32|) . v = 1. Let us use Eq. (|35|) to relate the asymptotic behavior to the 
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number of observation N and the analyzing scale index j. Since v = 1, we get that the 
asymptotic behavior of the scalogram S nj j depends on whether, as j,N — > oo, 

N < 2 2j or if 2 2j < N . 

Let us explain how these two regimes show up in the limit. The wavelet coefficients of Y can 
be expanded as follows 

where Wyjt , belong respectively to the chaos of order go and qo + 1. Then, 



The term [w\ q ^] 2 behaves as in the case G = H qo and is asymptotically Rosenblatt as 



= H?] 2 + [<r +1) ] 2 ) + (aw 



(<?o) w (go+l) 



3 



proved in IClausel et aD |2011bl ]. The term [W^ +1) ] 2 is asymptotically negligible as proved in 

Proposition 15.11 The term Wj^Wj^ , on the other hand, turns out to be asymptotically 
Gaussian. The asymptotic behavior of the scalogram then depends on whether the Rosenblatt 
term or the Gaussian term is leading. This depends on the limit of N/2 2 i . Hence, both limits 
stated in Theorem 13.11 may occur: 

• If 2~ 2jf N —7- 0, the term corresponding to [W^ '*] 2 is leading and the scalog ram S n ^ j 
of Y is asymptotically Rosenblatt. 

• If 2~ 2 i N —7- oo, the terms corresponding to VF^Vl^'? ^ are leading and the scalo- 
gram S Uj j of Y is asymptotically Gaussian. 



Example (iii) Suppose G = H qo + H qo +i + H qi with qo > 2, qo + 1 < 1/(1 — 2d) and 
Qi ~ (Qo + 1) > 2. Then / = {qo}, J = {(qo,qi), (qo + l,<7i)}- Observe that in this case, J 
is not involved in the limit of S nj j and the behavior of the scalogram is similar to that of 



Example (ii) Thus, the two limits of Theorem 13.11 may occur. 

3.2. The Hermite rank of G equals 1. Here we assume that 

q = 1, 3 < gi < 1/(1 - 2d) and £ = oo . (36) 

In particular, this condition implies d > 1/3, thus d £ (1/3, 1/2). By definition of Iq i n (|25p . 
the last condition in (j36[) means that there are no terms with consecutive indices in the 
Hermite expansion. Thus 

— -, H qi H 



G — c\H\ + — t H q ^ + —H q2 + 



where for any I E C, qe+i —qi>2. In this case the following critical index plays an important 
role : 

*' = 2 + W^J)- (37) 

It will also be useful to relate the number of available wavelet coefficients n = rij to 7J where 
v takes the following three values : 

(1 - 2d)( Ql - 1) I -2d gi-1 

Ul = l-(l-2d)( gl -l) ' V2 = M^Tj2 {qi ~ 1} ' " 3 = g7^3 • (38) 

As shown in the following lemma, the relations between ^1,1/3 and 1/3 depend on whether 
gi < g* or q x > q{ : 
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Lemma 3.1. 



• If qi < q\ then v\ < vi < U3 . 

• If qi > q\ then v-$<vi<v\. 

Proof. To prove Lemma 13. 11 observe that 



and 



(l-2d)(gi-l) 
l-{\-2d){ qi -\) 



< v 2 



(l-2d)(gi-l) 
2d - 1/2 



2d-- < l-(l-2d)(gi-l) 



qi < 1 + 



(l-2d)(gi-l) gi 
2d — 1/2 qi 



gi-3< 



2d-j 
1 - 2d 



9i < 3 + 



l-(2d-^) 
1 - 2d 

2d - A 



1 - 2d 



□ 



The next theorems indicate the limits in the various cases. We first consider the case where 
qi is lower than the critical index g|. 

Theorem 3.2. Suppose that Assumptions A hold with M > K+d. Suppose moreover that the 
Hermite expansion of G satisfies Ii36\) and assume that q\ < q*, where g* is defined in (31). 

Then three limits of the multidimensional scalogram S n j in U9\) . suitably normalized, are 
possible : 

(a) If Jj 3 "C rij then as j,nj — > 00 , 

■2 



-2(«(ffi)+J0q (£3 "gi 
n 3 1j "* J (gi-1)! 



(2tt)^- 1 [/*(0)] 91 L 9l _ 1 Z 2 , d (l 



fo/ <C nj <C 7J 3 £/ien as j, n.,- 



00 



(l-25(gi -l))/2 -(2tf(££i)+M0= (£) 47TCiC gi 



7 



(cj // nj <C 7J 1 i/ien as j, rij 



(91 " 1) 



00, 



1/2 -(2d+2i<r)c (C) 2 



'"3 ,3 



r, 



where T is defined as 
(/*(0)) 2 



7T 



I A + 2p7rr 2 ( K+d )[/i ii00 ^ i00 ](A + 2 P tt) 



pel 



dA, \<i,i'<m. (39) 



Remark 3.1. In case (a) the limit is a deterministic vector times the non-Gaussian Rosen- 
blatt random variable Z\^(X)- ^ n case [(b)] the limit is a deterministic vector times a Hermite 
random variable of order gi — 1>3 — 1 = 2, which can be represented by a multiple Wiener 
integral of order 3 or more (see Definition 12. ID . 

In the case where rij ~ ^o7j 3 as j — ?■ 00 for some Co > 0, the scalogram is asymptotically 
a linear combination of a Rosenblatt and a Hermite random variable. This is because it is 
the sum of two terms both converging in L 2 after normalization (see Section [6]) . On the other 
hand if rij ~ Corf 1 as j — > 00 for some Cq > 0, the situation is complicated. This is because 
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the scalogram is the sum of two terms of same order, one converging in I? to a Hermite 
random variable, the other converging only in law to a Gaussian random variable. 

Proof. This theorem is proved in Section 17.21 □ 

We now consider the case where q\ is greater than the critical exponent q\ . 

Theorem 3.3. Suppose that Assumptions A hold with M > K+d. Suppose moreover that the 
Hermite expansion of G satisfies 136\) and assume that qi>q\, where q\ is defined in (31). 

Then two limits of the multidimensional scalogram S n j in t!9\), suitably normalized, are 
possible : 

(a) If rij *C 7J 2 then as j,nj — > 00, 

nj /2 ^ d+2K %^ c?A/-(0,r), 

where T is as in Theorem \3.S\ (a). 

(b) If rf 2 <C rij then as j,rij — > 00, 

2 

n j Jj Z>n,j ~> _ {27T) 11 [f (U)\ L, ?1 _l Z 2t d\,l) ■ 

Remark 3.2. As in the case of Theorem I3.2[ the case where rij ~ Co7j 2 as j — ?■ 00, seems 
quite complicated to deal with. 

Proof. This theorem is proved in Section 17.21 □ 



Exam ple. We now illustrate Theorem 13. 21 and 13.31 Our setting is still that of lMoulines et al 
(20071 ] as above. 

The memory parameter d is assumed to belong to (3/8, 1/2). Consider the case where 

G = H\ + H Ql , 

with 3 < qi < 1/(1 — 2d). We will prove in the sequel that the wavelet coefficients of Y can 
be expanded as 

Wj* = w$ + wjf , 

where W-^ is Gaussian and belongs to the chaos of order qi. Then, 

The empirical mean of the terms [W^^] 2 behaves as in the Gaussian case and is asymptotically 

Gaussian. The empirical mean of the terms [Wj q ^] 2 behaves as in the case G = H gi with 

qi > 2 and is asymptotically Rosenblatt. Finally the empirical mean of the terms 2W^W^ q ^ 
belongs to the chaos of order qi — 1 > 2. The asymptotic behavior of the scalogram then 
depends on which of the three terms is leading. 

To see what happens, let iV be as before the number of observations and assume that 
7j = 2 3 . Let rij ~ N2~ 3 as j — > oo as in (|34l) . Distinguish two cases : qi < ql and qi > ql 
where ql is defined in ([37j) . 

If qi < ql, the three possibilities stated in Theorem 13.21 can occur : 

• if 2~ j( - ul+v >N — > as N, j — > 00, then the term corresponding to [W^] 2 is leading 
and the scalogram S nj j of the process {Y t }t^z is asymptotically Gaussian (case (c)). 
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if 2 -J ( l ' 1+1 ) N —7- oo and 2~^ V3+1 ^ 'N — > as N, j — > oo, then the term corresponding 
to SWjkWjk is leading and the scalogram S nj j of {Y t } belongs asymptotically to 
the chaos of order q± — 1 > 2 (case |(b)j ). 

if 2~i( u ' A+1 } N — ± oo as N, j — > oo, then the term corresponding to [Wj^] 2 with q\ > 3 



is leading and the scalogram S n -j of {Yt} is asymptotically Rosenblatt (case (a)) 



If we now assume that qi > q\ , we are in the setting of Theorem 13.31 and the term corre- 
onding to ^Wj^Wj^ is always negligible. Then only two different situations can occur : 

• if 2~i( U2+1 )N — >• as N,j — > oo, then the term corresponding to [W^] 2 is leading 



and the scalogram S n . j of {Yt} is asymptotically Gaussian (case (a)). 



— > oo as N, j — > oo, then the term corresponding to [Wj^ i ls lading 



and the scalogram S nj j of {Yt} is asymptotically Rosenblatt (case (b) ) 



4. The basic decomposition 
Our goal i s to investigate the a symptotic behavior of S Hj j as defined in f)19|) when j 



+oo. As in IClausel et al.l [201 lbl ]. our main tool will be the Wiener-Ito chaos expansion 



of S nj j which involves multiple stochastic integrals I q , q = 1,2,.... These are defined in 
Appendix IA1 In this case, the sit uation is more complex than in the case G — since 
as proved in IClausel et all |2011al ]. the wavelet coefficients W^jt, defined in <\18b . admit an 



expansion into Wiener chaos as follows : 



oo 

W ^ = £ 5 W S . ( 4 °) 

9=1 

where is a multiple integral of order q. Then, using the same convention as in (|20p . we 

have 

w?,* -£(»)■ + 2 EE B« • m 

g=l yH ' y q'=2 q=l H ' H ' 

where the convergence of the infinite sums hold in L 1 (il) sense. 

Each wj 9 ^ is a multiple integral of order q of some multidimensional kernel fj^ , that is 

wg = ?,(f$) ■ (42) 

Now, using the product formula for multiple stochastic integrals (|125p . one gets, as shown in 
Proposition 14. II that, for any (n, j) G N 2 , 

^ n— 1 

s—-E w ^- E i w io] 



'11,1 

n 



k=0 

oo / \ 2 9—1 / \ 2 



2(9,9,2?) 



E/Q ^ 

00 9'-l 9 / \ / /\ 

^EE^Ew^jQsW", (43) 



q'=2 q=l * * p=0 
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l(g,Q',P) ; c 



where, for all q,q' > 1 and < p < min(q, q'), S^'- is of the form 



o(q,q',p) _ f , c_.(9,g'>p)\ 

b n,j - I q+q'-2p{S n J ) ■ ! 1 ! ! 



We call q + q' — 2p the order of the summand S^ 9 For any n,j,q,q',p, the function 
gi7'' p) (0, £ = (6, • • • , W-2 P ) 6 R<?+4'-2p is defined for every p, q, q> as 

.. n— 1 

^'' p) (e) = ^E( f SM?) > ( 45 ) 



fc=0 

where the operation ® v is defined in (|126p for each entry. The expansion in Wiener chaos of 
S n j implies that 

„2 C (1,1,0) , v,(0) . v (i) . v,(2) . „(3) 



with 



Snj — c lS n ,j ' + S n,j + ^nj + S n,j + S nj > (46) 

S S = E t^E^{^ 2(2 ^ s -" p) ' (47) 



fc£,f>m y£ ' p=0 W VF/ 

-^(^^a^'')- <49) 

^^f>(:)(% + >M<: 



The sets £, / and J are defined in (f2Tj) . ([23]) and ([24"|) respectively and the index mo, 
defined in (j26|) . is such that g mo > 3. 

Let us comment on the decomposition (1461) , The sum E^- contains terms of the form S^' p ^ 
that is multiple integrals of order 2(q — p). Then this sum, after subtracting its expectation, 
has only summands of order 2,4,6,... in the Wiener chaos. 

The sum E^- contains multiple integrals of orders q + q' — 2p with q ^ 1, q' ^ 1, p < q A q' 

and \q — q'\ > 2. That means that all the summands in E,^- are of order greater than or equal 
to 2. 

(2) 

The sum E^ ■ contains multiple integrals of orders q + q' — 2p with q = 1, q > q mo > 3 and 

(2) 

p = or 1. All the summands in E^j are then of order greater than or equal to q mo — 1 > 2. 

The last sum E^- contains terms of the form S^j +1 ' p \ that is multiple integrals of order 
q + (q + 1) — 2p = 2q + 1 — 2p. When p = q, q+ 1 + q — 2q = 1, thus one can have components 
in the first Wiener chaos, that is Gaussian terms. 
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We will see that \ + S„ will converge to a non-Gaussian limit, more precisely to a 

random variable in the second Wiener chaos. The sum will also converge to a non- 
Gaussian limit, more precisely to a random variable in the Wiener chaos of order q mo — 1. 

(3) 

Finally will tend to a Gaussian limit. 

(2) 

Remark 4.1. It is the presence of j which creates the possibility of having as limit a 
multiple integral of order greater than 2. Thus, starting with a process 

G{X t ) = H l (X t ) + H qi (X t ) , 

with q\ > 4, then q mo = q\ and one may obtain as limit of the scalogram a Hermite process 
of order q\ — 1 > 3. 

Let us formalize the above decomposition of S nj - and give a more explicit expression for 
th e function gj 9 ' ? ' PJ in (145p . The next proposition is a generalization of Proposition 6.1 
of IClausel et aiT[2011bl ] . 

Proposition 4.1. For all j, {Wj^jfcg^ is a weakly stationary sequence. Moreover, for any 
{n,j) G N 2 ; S n j can be expressed as where the infinite sums converge in the L 1 (0) sense. 
The function g { n g f' p) (0, £ = (6, • • • , W~2 P ) G M 9+g '~ 2p ; m 0, equals 

E ( Z 9 '' P \0 = ^n(7,{6 + ---+W-2 P })xnS'" 2 nA/7(^l(-^)te)] (51) 

X «/ (£l H 1" £g-p>£q-p+l H 1" £,q+q'-2p) ■ 

Here f denotes the spectral density of the underlying Gaussian process X and 

i n— 1 -i inu 

DJu) = -Y e iku = . . . , (52) 
v 7 n ^ n(l-e m ) v 7 

denotes the normalized Dirichlet kernel. Finally, for £1,^2 £ ^ 3 i/p 7^ 0, 
«f(a,6) = ^ )p (ft^)) hf + - • - + A p +ei)nf } (A! + • • • + A p - 6) d p A , (53) 
and, if p = 0, 

^ P) (ei,6) = hf ) (ei)hf ) (6). (54) 

Notation. To simplify the notation, for any integer p and qi, ■ ■ ■ ,q p £ Z + we shall denote by 
the C qi+ - +q >> -»• C p function defined, for all y = (yi, . . . , y qi +... +qp ) G C« l+ - + * by 

(9i 91+92 <?iH hgp \ 

i=l 1=5!+! i=?iH — hgp-i+1 / 

Note that, for p = 1, one simply has S g (y) = y\ H 

With this notation, (|4"4"|) and ([53]) become respectively 

Si 9 / ' P) = W-2 P (A. o S 9+ ,'-2 P (7, x •) x [v7l(-.,.)] 0(<?+9 '- 2p) x fc« o E ? _ M ,_ p ) , 

(56) 

«f(6,6)= / /^(A) hf } (£ p (A) + 6)hf } (£ P (A) - 6) d p A , if p + 0, (57) 

J(-7T,7r)P 
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where o denotes the composition of functions, A = (Ai, • • • , A p ) and /® P (A) = f(Xi) • • • /(A p ) 
is written tensor product. 



Proof of Proposition \4-l[ For sake of simplicity we can assume that W is a vector of length 
m = 1 since the case m > 2 can be deduced by applying the case m = 1 to each entries. We 
must give an expansion in Wiener chaos of the one dimensional scalog ram S^i^j — in 
our setting. Using (fTUj) . (I4T1). (142]) and the pr oduct formula (I125P of Proposition lATl] we have 
as in Proposition 6.1 of Clausel et al. 2011bl ] 



^ = E f& E^w (!) (!) w-* • (^) 



where 

1 n— 1 

( M » _ I v f f (<?,) 

k=0 

By (M), 



= expoE^-O (h,™ o S ,(0) (f® q (Of 2 1^,(0 .^M 9 . (59) 

If g + g' - 2p ^ 0, let £ = (£i, • • • , £</+<;' -2p)- As in IClausel et all |2011bl ] using (fT2^D . we get 
that g^'j is a function with q + q 1 — 2p variables given by 



9nf' P \t) = i£ex P oE ?+g ,_ 2p (i& 7 ;£) x [V/W)] 8 ^-*® x *f ° £<™'-p(6 . 

fc=0 

The Dirichlet kernel D n appears when one computes the sum A X^fc=o ex P O ^g+g'-2p(i^7j0- 
This implies the formula (|5ip . 

In addition, the chaos of order zero appears in the expression (|58p of S^j in the ter ms with 
p = q = q' since I q+q i-2 P = lo- In this case, a similar argument as in IClausel et all |2011bl ] 
leads to 

^£7^£ E (KJ| 2 ) = £ W >< E d<| 2 ) =E(l^,o| 2 ) =E(S nJ ) , 

9=1 fc=l g=l W ' ; 

by (|30|) and ((HQ). Therefore, in the univariate case m = 1, S n j = S n j — K(S n j) can be 
expressed as stated in (fl3j) . The generalization to the case m > 2 is straightforward. □ 

We prove in Section [B] that 

c 2 

• The leading term of eJJ] + sJJ] is 7 ^ I ^(27r)9o*- 1 S^.' 9o *' <7o *" 1) (see Propositions O 
and !6.2p where 

^° | g otherwise. 

Note that sf, 9 °-' 9o ' 9 ° ^ always is in the 2nd Wiener chaos. 

• The leading term of E^. is 2ir ? 1 9m °,, S^ 1 ' 9 " 10 ' 1 ^ (see Propositions 16.31 and I6.4|) . which 
is in the (q mo — l)-th Wiener chaos. 
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• The leading term of is 2 ^^q11 Qtt^o S K^o +1 ^o) ( gee Propositions [63] and 
16.61). which is Gaussian. 



Hence, for the two classes of functions considered in Sections 13.11 and 13.21 we have to com- 
pare at most four terms : S^ 1 ' ^, S^-' 90 ' 9 '* 1 \ S^j' qt ° +1 ' P , S^ 9 ™ ' 1 ^ which are respectively 
asymptotically Gaussian, Rosenblatt, Gaussian and in the chaos of order q mo — 1. 

Our three theorems are based on the study of the asymptotic behavior of each sum (see 
Section [6] below) . We first establish some preliminary results. 

5. Preliminary results 



5.1. L 2 bounds. To identify the leading terms, using the same approach than in lClausel et al 



|2011bl | we will give an upper bound for the L norm of the multidimensional terms Sir'? , 
q,q',p defined in (PHI) and (|56l) . Here, the main difficulty is that unlike the case where G = H qo , 
w e have to deal with an infinity of terms. We have also to obtain more precise bounds than 
in IClausel et al. 2011b ], In the following, for any random vector Z, the L 2 (VL) norm of Z is 
denoted by 

||Z|| 2 = (E[|Z| 2 ]) 1/2 . (61) 

(Recall that |Z| denotes the Euclidean norm of Z.) Our goal in this section is to specify how 

||Sif',' ' H2 depends on q,q' and p. The difficulty is that the sum Syj j contains long— memory 
and short-memory terms having then different normalization factors. To recover all the cases, 
we shall use not only 5+(q) and 5(q) defined in (|12p but also 

S-(q) = max(-%),0) , g>0, (62) 

so that 5 = 5^—6- and <5+,(5_ are nonnegative. In particular, 5(0) = <5+(0) = 1/2 and 

<5_(0) =0. 

As in Clausel et al. 2011b ], the expression (|56p of S,^ 9 ! 9 involves the kernel defined 
in (157p and we have to distinguish the two cases p / and p = 0. The following notations 
will be used in the sequel. For any s G Z + and d G (0, 1/2), set 



A s (a) = H(a i [) 1 - 2d , Va = (ai, • • • , a s ) G N s . (63) 

For any q, q',p > 0, set 



a(q,q>,p) = r in ( 1 - 5+iq - p) - 6+iq '- p) >^ if ^°' (64) 
[2 if p = , 

fi(q,p) = max \6+(p) + S + (q - p) - ^,oJ . (65) 

Notice that for any q > 0, (3(q,0) = 5+(q). Define the function e on Z + as 

e( P ) = l° ifforan y se {V" ) i'}, Ki- 2 ^ 1 . sqqs 

I 1 if for some s G {1, • • • ,p}, s(l — 2d) = 1 . 

The index K is defined in ([3]) and the index M is defined in (|129p . and, as noted in AppendixlBl 
the filter hj(t) has null moments of order 0, 1, . . . ,M — 1. 

Proposition 5.1. Suppose that Assumptions A hold. Let e be the function defined in (E 
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(i) There exists some C > whose value depends only on d and f* such that for any 
n, j > 2 and any < p < q < q' 



Si 9 /' p) || 2 < C 3 ^A 2 (q-p,p) 1/2 A2(q' -p,p) 1/2 n- a ^'' P Hlogn) 



e{q+q'-2p) 



X7 



2K+l3(q,p)+f5(q',p) 



J 



(logjj 



|3e(g') 



(67) 



(it) In addition, assume that M > K + max(<5+(g), 5+(q')). Then there exists some C > 1 



whose values depend only on d and f* such that for any q,q ,n, j 



(68) 



Proof. As above, we can take m = 1 without loss of generality. In what follows, C, C\,- ■■ 
are positive constants that may change from line to line. The following function e' defined on 
M+ is used in the sequel, 



e'(a) = 1 {1} (o) 



e'(a 



1 if a = 1 , 
s'(a) = otherwise. 



(69) 



We first observe that if p = 0, the bound ([67]) is a consequence of ([68]) . This follows by 
observing that A2^q, 0) = Ai(g), a(q,q',Q) = 1/2, 0) = and that the log exponents 

in (|67|) are all larger than that in (|68p . Hence, to prove the result, we show that (i) holds for 
p > and then prove (ii) , successively. 

Proof of (i) for p > 0. Set r = q — p and r' = q' — p. The starting point of the proof is the 
integral expression of given by (|56p . Thereafter we follow the same approach as in the 

proof of Proposition 7.1 of Clausel et al. 2011bj ]. using Lemma l8lfl to bound the kernel Kj 

involved in the integral expression of S { ^f' p) instead of Lemma 10.1 of lClausel et"afl [201 lbj |. 
replacing 2r, (r,r), 5{p) with r + r', (r,r') } S+(p) and adding if necessary a logarithmic 
correction. 

We obtain the following inequality, similar to (7.2) and (7.3) in lClausel et al.l 201 lbl ]. 



E 



S 



(<?><?',?>) 



where, for any j, n, 



7j7rr 



IP 



IP 



< C 1 2p (p!) 2 ( 1 - M ) 7 7 2+25(r)+25(r ' )+45+{p) 7f (log7i) 2£(p) /n,i , (70) 



r> (1+n | {t4l+ ^ } | )2 (l +7j 



25+ (p) 



1 + 7; 



l 7j ' 



28+(p) 



and where J r ,-yjir{u\\ 2dl r ) and ^',7,^(^1; 2dl r i) are defined in Lemma ED 

We now use inequality (|105p of Lem ma IHTTl succes s ively w ith p = r, a = 7 j7r, s\ = ui and 



p = r' , a = 7j7r, s\ = v\. As in (7.4) in lClausel etaD [2011b| |. we get that 
Inj < C 2 r+r '(r!r / !) 1 - M 7 f- (r)+2 ^ (r ' ) (log7i) £(r)+£(r ' )+2£(p) 



(— 7rr,7rr) ( ■y J - )^(— ivr' ,ivr') \ j. 



) ( ^ ) i«i r 2<5+ (r) 1 «i r 25+ (r ' ) d«i dui 



(l + n|{ui+«i}|)^l + 7 j {f} 
where C2 > 1 denotes a constant depending only on d. 



25+ (p) 



1 + 7, OJ} 



2<5+(p) 
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The next step relies on the inequality \{x}\ < \x\ on M. and on the 2-7r~periodicity of x i— > {x}. 
We then get that 



with 



Inj < C 3 r + r \r!r / !) 1 - 2 S7 25 - (r) - 25 - (r ' ) (log 7j £(r)+£(r ' )+2£W in,i , 

| Ul |-2*+(r)| t;i |-M + (r') dui dvi 



l(- 7j ^)2 (l + n\{u l + v 1 }\)*(l + \u 1 \)M+<P)0. + \vi\) 2S +W ' 

which corresponds to (7.5) in lClausel et al. |2011bl |. 

The bound of I n j is obtained using the decomposition J n j = A + 2B with 

| ni |-25+(r)U 1 |-2 1 5 + (r') dui dvi 



A 



L m (1 



and 



where 



s=l 



+ n \{ U1 + wi}|) 2 (l + |ui|) 2 Mp)(1 + \ Vi \) 2S +(p) 



A- 



« (l + n|{m + fi}|) 2 (l + \u 1 \) 25 +(p){1 + \v 1 \) 2S +(p) ' 



G (-7j7r,7 : ,-7r) 2 , |ni + di - 2irs\ < ir} 



with s G {~ 7jj • • • ,7j}' Th i s decom position is similar to the one used in the proof of Propo- 
sition 7.1 in Clausel et al. [2011b ] and is obtained by partitioning (— 7j7r, jjir) 2 using the 



domains A^ . 

We now bound separately A and B, as in I Clausel et al. 2011b ]. We get that there exists 
some C > such that if 28+ (r) + 2(5+ (r') > 1, 



and if 25+{r) + 25+(r') < 1, 



A < Cn' 1 (log n 



'(25 + (r)+25 + (r')) max(l-2<5+ (r)-25+(r-')-45+ (p),0) 
/ 1 



(log 7j 



.^(2«+(r)+2« + (r')-HM+(p)) 



where the function e' has been defined in (16911. Further for some C > 



5 < ^-iy^ 1 - 2 ^)- 2 ^)' ^^ 

As in I Clausel et al. I [201 lbl ], we deduce that there exists some C > depending only on 
5 + (r ) , <5 + (r') , 5+ (p) , d such that 

J nj < (7 n -mm(2(l-5 + (r)-5 + (r')),l)( lo g n )e'(25+(r)+25 + (r')) 

max(l-2(5 + (r)+5 + (p)),0)+max(l-2(5 + (r')+<5+(p)),0) n ^'(25+ (r)+25+ (p))+e'(25+ (r')+25+ (p)) 

x 7j - l^iog 7j ; 

Observe now that for any fixed ci, there exists only a finite number of possible values for 
5+(r),5+(r'),5+(p) and then a finite number of possible values for C. Then, provided we 
replace C by its maximum possible value, we can assume that C does not depend on <5+(r), 
5+(r% 5 + (p). 
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The bound on I n j and ([7D]) yields 

\\S { ™'' P) \\ 2 <C q -^A 2 (q -p,p) 1 / 2 A 2 (q' - p ,p)V2^K n ->*Hl-S+to-p)-S + W-l>)A/2) 
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x 7 



max(l/2-5 + (p)-<5 + (g-p),0)+max(l/2-5 + (p)-5 + ( g '-p),0) 
-1+5+ (q-p)+S+(q' -p)+28 + (p) 



X 7 , 

X (l ogn )e'(2<5+M+2<5+(r')) 

(log7 i ) [e ' (25 + (r)+25+(p))+e ' (25+(r ' )+25+(p))+E(r)+£(r ' )+2e(p)]/2 , 

Inequality (f67"|) corresponds to this bound with exponents of 7 j, logn and log 7 j simplified as 
follows. 

The exponent of 7 j is obtained by observing that — 1 + d+(q — p) + 5+(q' — p) + 25 + {p) = 
(—l/2 + 5 + (q — p) + 5 + {p)) + {—\/2 + 5 + {q' — p) + 5 + {p)) and using max(— a, 0) + a = max(a,0) 
with a = — 1/2 + <5+(g — p) + 5+(p) and a = — 1/2 + 5 + {q' — p) + 5 + (p) successively. 

The log exponents are obtained by observing that, since r < r', e(r) + e(r') + 2e(p) < 
2(e{r') + e(p)) < 4e(r' Vp). In addition e'(25+(m) + 2<5+(m')) = iff m + w! / 1/(1 - 2d) 
and equals 1 otherwise. Thus e'(25+(m) + 25+ (m')) < e(m + m') and we get 

e'(25 + (r) + 25+(p)) < e(q), and £'(25+(r') + 25+(p)) < e(q') . 

Finally, since e is non-decreasing and q < q' ', r' V p < q' ', 

e'(25+(r) + 25 + {p)) + e'{25 + {r') + 25+ (p)) + 4e(r' Vp) < e(g) + e(g') + 4e(g') < 6e(g') . 

Proof of (ii), Here, p = and t hus = hj^ 2 . The same approach as in the proof of 
Proposition 7.2 in IClausel et alJ 201 lbl ] leads to the following inequality which corresponds 
to (7.12) in IClausel et al.l [201 lbl ] : there exists C > not depending on n,j,q,q' such that 

n\s { n q f' 0) \ 2 } < c 7 ; {q+q,){1 - 2d \f 2K+1) i nj = Cl fww)+^) In . i (71) 

where 



q-yj-K 



-Tiji 



q 7j7 



g(u,v)J q>ljn (u;d,--- , d) J q > , lj7r {v, d, ■ ■ ■ ,d)dudvi , 



(72) 



with J m ,a defined as in Lemma 18. II and with g(u,v) defined for all (u,v) G R 2 by, 

g(u,v) = (1 + \n{u + ,}|)- 2 \^ U l^\ m ~ K) • ^ {vh ^ M ~Z , . 

As in the case p / 0, we can use the bound (|105p of J m>a and the inequality \{u}\ < \u\. 
We get that 

I n ,j < C^^qlq'l) 1 - 2 ^/^- 5 -^ 



97.7 7 



J 

q'-yj-K 



7i{£} 



2M-2K-26 + (q) 



7;{f} 



2M~2K-25+(q') 



dudv 



2{M+a) 



2(M+a) 



As in the proof of Proposition 7.2 of IClausel et al. [201 lbl ] . we then obtain that 

Inj < Cl +q \q\q'^- 2d n-\ 2{S - {q)+S - m . (73) 
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The conclusion follows from (|7ip and (|73p . □ 

5.2. Asymptotic behavior of the leading terms. We now investigate the exact asymp- 
totic behavior of the terms that will turn out to be leading in the sum (143|) . 

Let us first suppose that the bounds in Proposition 15.11 are sharp enough to determine 
which terms are leading. Since jj — > oo and n = rij — > oo, those for which the bounds 
have the largest exponents f3(q,p) and /3(q',p) and the lowest exponent a(q,q',p) are more 
likely to dominate, in particular, if 5 + (p), 5 + (q — p), 5 + (q' — p),5 + (p) + S + (q — p) — 1/2, 
5+(p) + d~ + (q' — p) — 1/2 and 1/2 — (1 — 5 + (q — p) — 5 + (q' — p)) are all positive. Using (fT5j) . 
if p > 0, this happens for < p, q — p, q' — p, q, q' , q — p + q' — p < 1/(1 — 2d), that is (taking 
q < q' without loss of generality) , 

0<p<q<q' <l/(l-2d) and < q + q' - 2p < 1/(1 - 2d) . (74) 

In particular, for such a triplet {p,q,q'), we have e(q') = e(q + q' — 2p) = so that bounds 
in (I67p and (|68p involving logarithms will not appear in these terms. We shall check afterwards 
(in Section [6|) that indeed, in all the cases we consider, either such a term is leading in the 

sum (|4"3j) . or the leading term is S^ 1 '* (q = q' = 1 and p = 0). The bounds established in 
Proposition 15. II will be sharp enough for this goal. 

This is why, in the following, we shall only determine the asymptotic behaviors of S^7 1 J 

and of S^.'j' 3 ^ under Condition (f?4"|) . when j,rij — > oo. 

Proposition 5.2. Suppose that Assumptions A hold with M > K + 5(1) = K + d and that 
7,- is even for all j. Let (nj) be any diverging sequence of integers. Then as j —> oo, 

nf^^^M^T), (75) 

where V is defined by L39\) . 



Proof. This is a direct application of Theorem 3.1 case (a) in lClausel et al.1 [2011bl |. □ 

We now consider the case where Condition (|74p is satisfied. 

Proposition 5.3. Let q, q' and p be non-negative integers such that ( |7^| ) holds. Assume that 
Assumptions A hold with M > K and let (nj) be any diverging sequence of integers. Then, 
as j — > oo, 

( ni7i )l^(a-p)^( 9 '-p) 7 7 2 (^+ 5 W)s^^ ( 4 } [r(0)]^')/ 2 L p Z q+q ^ 2p4 (l) , (76) 
where Z q+q >^2p,d * s the Hermite process defined in \2$fy and L p is defined in \33\) . 



Proof. The proof follows the same line as the proof of Proposition 8.1 in lClausel et al.l 2011bl |. 
Therefore we only explain how to adapt this proof to our setting. Set r = q+q'—2p. Using (|56p 
and that, for all g G L 2 (W), 

1(g) I (n l3 rl 2 l(g(./(n l3 ))) , 

we have 

S?/' P) = (ni 3 )- T/2 % (D n o ^ q+q <-2 P (-/n) x [l ( _ 7j w) r(-/n) x £,) , (77) 
where, for all ^Gi r , 

fi(«7i0 = ^ r \o x *f o ^ q - P , q >- P (o ■ 
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The rest of the proof consists in proving the 1? convergence of the ltd integral i n (1771) . 
adequately normalized. This is done in the proof of Proposition 8.1 in lClausel et al. 2011bl | 
with q — p = q'— p = 1 (hence r = 2). The same proof applies in our setting but results in 
a multiple integral of order r with r > 2. In particular, if r > 2 the asymptotic limit is not 
Rosenblatt but an r-order Hermite process. □ 



6. Leading terms 



Recall the decomposition (|^6|) of S n j using sums aim of 

iection is to identify the leading terms of the three follo\ 
^njj un der the conditions specified in Sections 13.11 and I3.2L 



this section is to identify the leading terms of the three following sums : E^ • + E^/ -, E^/ ., 

11 j, 3 ll j>3 li j'3 



6.1. Leading term of E^ ■ + E„ „•. Recall that the two sums „-,e£^ • are defined in 

_ n j ,3 < L j ,3 li j i J n j i J 

equations (|47|) . (j48l) and that q$, defined in ([60]) . equals q\ if qo = 1 and equals go otherwise. 
Therefore go > 2. We shall prove that, if < 1/(1 — 2d), the main term in Ej^. ■ + E^ ■ is 

( i ^ ITT (27r)«o*- 1 s[ l 9 °' <?0 *'' ?0 *" 1) , and has rate n l-2d 7 ~2(%5)+JO > Xhe following proposition is used 
to show that the remainder terms are negligible. 

Proposition 6.1. Suppose that Assumptions A hold with M > K + 5(q$) and that 

q* < 1/(1 - 2d) , 

where q^ is defined in Let (nj) be a diverging sequence. Then, when j —> oo, 



11 : 



2 / *\ 2 



7r wor ; >: t^pv: w p n s n7,i iia -> , (78) 



1-^-^,,,-,1+a, [ V _^_^ p! W (27r)P||S^f' p) || 2 ^0, (79) 

taC, qi >q* wv p=0 / 



,1-M -2(«(9S)+#) 
l i 1 3 




(80) 



Proof. We first note that, since > 2 by definition and q$ < 1/(1 — 2d) by assumption, we 
have 1/4 < d < 1/2. 

As in the proof of Proposition 15.11 we prove the result in the case where m = 1 without 
loss of generality. We thus use non bold faced symbols. 

We first prove (|78|) . Since there is a finite number of terms in the sum appearing on the left- 
hand side of ([751) . it is sufficient to show that each term converges to 0. Let pG{0,...,gJ-2}. 
We apply Proposition 15.11 with q = q' = q%. Since q$ < 1/(1 — 2d) and thus e(q$) = 0, 
Inequality (|6T|) reads 



7 -2(%*)+A) ||5 ^, (?S , P ) ||2 < c ,SA 2 ( g 5 -p,p) n - a ^p) ^ ogn fm-p)) 7 f^o^-^o)) 

By ([109]) and lfIT5|t in Lemma E2\ we have a(qfi,q$,p) > min(2(l - 2d), 1/2) and P(q* a ,p) < 
*+(?o) = %o)- Hence > 

7 T mq ° )+K) \\S^' q °' p) \\ 2 < C c >oA 2 (q* -p,p)n-^ 1 - 2d ^ (lognf 2 ^-^ . 
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Since d G (1/4,1/2), we have n -min(2(l-2d),l/2) (l og n)< 2 ^o-v)) = {n 2d - 1 ). Thus Inequal- 
ity £78]) holds. 

We now prove (|79|) . We apply ([67|) . in the cases q = q' = qg, < p < qg — 2 and q = q' = qg, 
p = qt — 1, successively. Then for some C > 0, for any I G C such that ^ > q$ and for any 

< p < % — 2 one has, 

7 -2(% *)+if) ||5 fe, %1 p) ||2 < C ^K 2 {qt-p,p)n-^^ (log^^^^-^^f^^^-'^^aogT,) 3 • 

On the other hand, since d > 1/4, e(2) = 0. Thus in the case where p = qg — 1, the exponent of 
logn vanishes. Moreover, in this case, by (jllOf) in Lemma f8.2] a(qg, qe,p) = min(l — 2d, 1/2) = 

1 — 2d. In the alternative case p < qi — 1, we use (|109j) in Lemma 18.21 which gives 

( logn ^(2(*-P)) < n -mm(2(l-2d),l/2) ^ g n j 6 (2(«-p)) < n M-l _ 

for n large enough, since 2(1 — 2d) > 1 — 2d and 1/2 > 1 — 2d. Hence in all the cases, the 
terms in n can be bounded by C'n 2 ^ 1 . As for the terms in jj, we use that, by (I112p in 
Lemma 1831 f3(qg,p) < £+(%) < 5+(go + 1)> since > t/g + 1 and 5+ is non-increasing. Hence 
we get that 

7 -2(% S ) + ^) ||5 ^, 9£lP)||2 < ^ Aa( ^ _ P)P)n 2 d -i ^w+iHW))^^ > 

where C > may have changed from the previous line. Hence, summing over I G C such that 
% > q$ and all p E {0, ...,<# — 1}, we get that 

eEt&Kp) wings-"!. 

< if M9S+i) - %s)) (io g7j ) 3 ^ ^ (a^co* e p [ ( ) A ^ • (si) 

By definition of A2 in (|63p . for any integer q and any < p < q, 

p\(^\\ 2 (q-p : p) = (p!)- JW (g!) 2 ((9-p)!)- 1 - 2d < (gOV)" 2 * • (82) 

Hence the sum over p in the right-hand side of ([81]) is bounded by (qgl) 2 up to a multiplicative 
constant. In turn, the sum over I > 1 is bounded, up to a multiplicative constant, by 
Ylt=i c qe (2itC) qi , which is finite by Condition ([29]) . Finally we observe that, since 6(q^ ) > 
<5(<7o + 1) and <5((?o) > 0, we have 

7 ^ + ( 9 o+l)-%o)) (log7 . )3 ^ 0) 

as 7j — >• 00. Hence Inequality ([79]) holds. 

We finally prove that (f80|) holds. Inequality (f67|) for (^1,^2) £ ^ with q = qe 1 ,q' = qt 2 and 
P < % a implies that 

-2(6(qZ)+K),, a (q ei ,qt 2 ,p)u ^1+^2 1/2 1/2 

7j ll^n,/ lb < 2 A 2(% -P,P) 7 A 2 (% 2 -P,P) 1 

x „-afe 1I?%)P ) (logn)£ (^ 1 +* 2 -2p) ^(^ 1 ^)+/3(^, P )-25(< ? 5) (log7j) 3 (g3) 

We first bound the terms that depend on n. First suppose that p = qg 1 and qg 2 = qg 1 + 2. 
In this case, the exponent of logn vanishes, since e(2) = for d > 1/4, and by (|l(J7j) in 
Lemma [521 the exponent of n a(qg 1 , qi 2 ,p) > 1 — 2d. Hence, in this case, the terms in n are 
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bounded by n 2d_1 . Otherwise, if p < or qg 2 > qi x + 2, we observe that for (£1,^2) G <A we 
have p < qe 2 — 3 and hence, by definition of a in fj64f) and since 5+ is non-increasing, 

p) > l/2-<5 + (3)=min(l/2,l/2-(3d-l)) > 1 - 2d , 

since 1/4 < d < 1/2. Whatever the exponent of logn, we again obtain that the terms in n 
are bounded by n 2d_1 , up to a multiplicative constant: 

sup n- a ^i^2.P)(logn) £fe i + ^- 2p) = C>(n M - 1 ) . (84) 

We now bound the terms that depend on jj in f|83f) . By (I112p in Lemma 18. 2[ we have 
(3(q,p) < 6+(q) for < p < q. Thus f3(q h ,p) + /3(%,p) - 5(q* ) < 6+(q h ) + 5 + (qe 2 ) - 2<%*). 
Since <5 is non-increasing, q^ > q$ and qi 2 > q^ + 2 we deduce that 5 + (q£ 1 ) < S+iq^) = S(q^) 
and 5+(qe 2 ) < 5+ (^ + 2) < £(<7o)- Hence the exponent of jj is bounded by a negative constant 
and 

sup 7» (l°g7j) ~~ ^ as j — )• oo . (85) 

(il,i2)eJ,0<p<q ei 

In view of ([55]). ([51]) and ([55]). the proof of (JgDJ) follows from the bound 

£ ^yC^^!(^]]W[A 2 fe-„)]^<oo, (86) 

(£ l/2)eJ «1- p=0 i=l 

which we now prove. By the Cauchv-Schwartz inequality, we have 

1 /9 

Qix 2 2 / 2 ^ 

n ( % W % -^ i/2 ^ n e(po 2 (2ttf ( % ) a 2 ( % - P ,p) 

We may replace % 1 by q^ in the previous upper bound, since qi x < q£ 2 . Likewise the sum 
over (ji,j2) £ J in ([86]) is bounded by the sum over (ji, j'2) £ N 2 . Hence ([86]) follows from 



, V2 

2 



M(qe-p,p)j < 00 , 

which is a straightforward consequence of (|82[) and the fact that, by Condition (|29[) . 

+00 

E|c % l(2vrC)*/ 2 <oo . 
t=i 

This concludes the proof. □ 

We now focus on the leading term of the sum X„ ■ + .. 
Proposition 6.2. Under the same assumptions as Proposition \6.1[ we have, as j —> 00, 



2 



V^ /,(0) '" L ' 5 - 1 



Z2,d(l) ■ (87) 
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Proof. We apply Proposition 15.31 with q = q' = and p = q$ — 1. Since 

25(1) + 26(q* - 1) - 1 = 2%5) , 

we get that 



„l-2d„-2(%g)+^) ^g ( 27r ) 9 ° 1 g (qS,qg,gS-l) (£) 



A (27t)%*- 1 



^2,d(l) • (88) 



The left-hand side in ([88]) corresponds to the term q% = q$ and p = q^ — 1 of in 
The terms of Xl[ l ' > ■ with q£ = <7q and p < — 1 are gathered in the left-hand side of ([75]) . 



The terms of j with % > t/g are gathered in the left-hand side of (|79p. Finally the left- 



hand side of (j80j) corresponds to - j m (SH) - Hence, by Proposition 16.14 all these terms are 



,(o) 

; wiLii m ^ yo e 

negligible and (f87l) holds. □ 

(2) 

6.2. Leading term of ^ -. In this section, we investigate the asymptotical behavior of the 
sum E„ j defined in (|49l) . We shall prove that, if q mo < 1/(1 — 2c£), the leading term of this 
sum is 2tt ^ g ^ ! S^ ?m °' 1 ^ and has rate ^ 2<5( -' ?mo 1 ))/ 2 ^%mo)+ c( + 2 - ft ^ rp Q gn( j wg 
show that the remainder terms are negligible. 

Proposition 6.3. Assume that Assumptions A hold with M > K + d and that 

q mo < 1/(1 - 2d) , 

where q mo is defined by [26\). 

Let (rij) be a diverging sequence. Then, as j — > oo, 



l/2-5(q mo -l) nr d{q mo )-d-2K [ C\C qi || C (l,gi,0) 



E ^HS^IIa Wo, (89) 



^>mo 



n l/2-S [qmo -l)^ S{qmo y d -2K 2n ^ ^1% || S (l.g.l) ||2 (90) 

\ e>m ^ qe j 

Proof. Observe that <5+(l) = d. We apply (|68j) in Proposition 15.11 with q = 1 and g' = qi. 
Thus there exists some C > such that for any £ > uiq 

7 - a ^o)-^|| 5 a^ ,o)|| a < cr^( %! )i/2-d n -i/^+(*)-*(*»o)(iog 7:7 .^(«) . (91) 

Since by assumption q mo < 1/(1 — 2d), we have e(q mo ) = and 5 + (g mo ) = 5{q mo )- Thus the 
terms involving jj vanish in the right-hand side of (|91|) . If i > mo, we have < 5(q mo ) 

and these terms are o(l) as j — )• oo. Hence, for j large enough, and for any t > mo, 

Using that 5(q mo — 1) > 6(q mo ) > 0, and that, by Condition 



we obtain the limit (|89|) . 
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We now show that ([90]) holds. Applying ([67]) with q = 1, q' = qi and p = 1, we get that 
there exists some C > such that for any ^ > mo, 

lisgf'^lla < C q -^{( qi - l)!} 1 /2^ n -(^^) (logn)£ fe-i)^(i,i)+/3fe,i)+2^ (log7 . ) 3 (Q2) 

The definition of a and /3 by Equations ()64p and (I65p . implies that 

a(l, ft, 1) = 1/2 - 5+(q e - 1), /3(1, 1) = d, p(q e , 1) = max(d + 6 + {q e - 1) - 1/2, 0) . 
Since £ > mo, one has 8 + (qt — 1) < 5 + (q mo+ i — 1). Thus 

-a(l,qt,l) n \e(ae-l) ^ 1/2-5+ (<?m n + i -1) , / l/2-5(g m , 

n- y '(iogrij) w ' <n- \ogrij = o[n- ° I . 

Observe now that for £ > mo, we have qi — 1 > q^o and thus 

y3(l,l)+/3fea)+21^ log7 ^3 < ^d+2K+m^(d+5 + (q mo )-l/2,0)^^ 3 = q f d+2K+S(q mo )^ 

Now, using the last two displayed equations, ([9"2"]) and Condition (|2"9"j) . we obtain the limit ([9"U]) . 
which concludes the proof. □ 

(2) 

We now deduce the asymptotic behavior of ^. 
Proposition 6.4. Under the same assumptions as Proposition I fi. 31 we /iaue as j — >• oo 

-W^)^-«0 S W (4) i^^[ / *(0)](-, + i)/ 2Ll Vo _ M (l) , (93) 



where Li is defined in k°3'3\) and Z q _\^ is the Hermite process defined in 

Proof. We apply Proposition 15.31 with q = 1, q' = q mo and p = 1. For these values, since 
<7m < 1/(1 ~~ 2d), Condition ([74"]) is satisfied. The exponents of n and 7j in the left-hand side 
of ([76]) respectively read 



l-S(q-p)- %' -p) = l- 5(0) - 5(q mo - 1) = 1/2 - 5(q mo - 1) 

and 

l-S(q-p)- 5{q' -p)-2K- 25(p) = -5(q mo ) - d - 2K . 
Hence we get that 



n . -^o-l))/2 7 -(5(^o)+^^) s a i , mo , 1 ) (£) [rm ( qmo+ l)/2 Li Vo _ m(1 ) . (94) 



(1-25( 3t , 

Finally we observe that this term corresponds to the second term of the summand in (|49p 
with index I = q mo , up to the multiplicative constant 47rcic (?mo /(g mo — 1)!. All the other terms 
are negligible by Proposition 16.31 Thus the limit (|93[) holds. □ 

(3) (3) 

6.3. Leading term of j- In this section we investigate the asymptotic behavior of j 
defined in ([50]) . We first bound the sum over indices £ = £0 and p / qi and the one over 
indices £ > £0 and p E {0, . . . , q{\. The two sums will turn out to be negligible. 

Proposition 6.5. Assume that Assumptions A hold with M > K + 5(qe ) and 

% + 1< 1/(1 - 2d) . (95) 
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Let (rij) be a diverging sequence. Then, as j —> oo, 

y- c ^o fgyH (2 _ )P f qe Q \ ( qe + 1\ ( 9v%0 +i, P ) 
^ % ! (% + !)! VP / V P / 



2 




Proof. Observe that, since qi > 1, the assumption g^ + 1 < 1/(1 — 2d) implies that d G 
(1/4,1/2). 

We first prove Inequality (|96p . Since there is only a finite number of terms in the left 
hand side of Inequality f)96h . we only have to prove that each term tends to 0. We apply 
Proposition 15.11 with q = qe , q' = qt Q + 1 and p < qi Q — 1. For these values of q, q' and 
p, under Condition (|95l) . we have e(q') = 0, and by (jlllj) and (|112l) . we have a(q,q',p) > 
min(3(l/2 - d), 1/2), < 8+(q lo ) = 5(qe ) and P(q',p) < 5 + (q, + 1) = % 4 + 1). Thus 

Equation (fBTl) yields 

n (l-2d)/2 7 -(%£ )+ 5 feo+ 1 )+ 2 ^)|| 5 ^V^O+ 1 'P)|| 2 = Q ^ n -min(l-2d,d) log ( n )^ 

Since d G (1/4, 1/2), we obtain (l96|) . 

We now prove (|97|). We apply Proposition 15. II with q = qp, q' = qg + 1 and p <qi for some 
£ £ I \ {£o}- In this case Inequality ([67|) reads 

7," IKj ||2 

< C^ + 2A 2 (%-p,p)U 2 (^ + l -p,p)ln- a («'« +1 'P) log(n) £ ^ +1 - 2 ^ 

x 7j Og7j) • (98) 

We observe that for n large enough, 

n -afe, % +l >P ) log ( n )«(29H-l-2p) < „-(l-2d)/2 _ (gg) 

Indeed, on the one hand, if p = q#, then e(2q£ + 1 — 2p) = e(q£ + + 1 — 2%) = e(l) = and 
a(qe,qe + l,qi) > (1 - 2d)/2 (fTUgj) . On the other hand, if p < q e , since d > 1/4, (fTD7|) implies 
that «((/£, + l,p) > 1 — 2d. 

In addition, by (|112|) one has for any p < q^, f3(qe,p) < <5+(%). Thus, for any i > £ and 
any p < q e , 

W(q e ,p)-5(q i0 )) + (P(qi + hp)S(qi +-L)) n x 3 ^ (<5+(%)-<5(% ))+(<5+(%+l)-5(g A) +l)) , 3 

lj (log7j) <7j (log 7i) 

< ^ + fe 0+1 )-5( % ))+(5 + ( %+1 +l)-5fe +l)) (log ^ )3 = o(1) _ (10Q) 

As in the proof of Proposition 16. 11 the Cauchy-Schwartz inequality and Condition (|29p imply 
that 

e (2i T?w" +1 Ep'( 9 ') ( 5 ' + > 2 (*- P ,P)*A 2fe+ i- P ,p)i < oo . 

q t \{q t + l)\ ^ \Pj\ P ) 
Applying this, O, & and ([TDD]> . we obtain $7$. □ 
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The following result can now be established. 
Proposition 6.6. Under the same assumptions as Proposition \6.5l we have as j — > oo 

n (l-2d)/2 7 -( 5 K + l)^fe ) + 2K) E (3^ (£) 2{27r)qeo [/» ( )«o+V2 Lgfo]ZM(1) 



Proof. We apply Proposition 15.31 with q = q' — I = qi and p = qe - Indeed we have, under 
Condition (|95j) . < g = qi < q 1 = qi + 1 < 1/(1 — 2d) and q + q' — 2p = qi + q£ + 1 — 2g£ = 
1 < 1/(1 - 2d). Thus Condition ((711) holds. We obtain that, as j -> oo, 

fl-acQ/2 -(5fe )+5K + l) + 2ir) 5 fe ,«„+!,«„) (£) [r(0) ^+i/2 Lft . (101) 



Using this limit, Proposition 16 . 51 and the definition of m dSOJ) » we conclude the proof. □ 

7. Proofs of Theorems EHJ GL2] and ET31 

7.1. Proof of Theorem 13.11 In the setting of Theorem 13.11 one has go > 2 and thus ci = 
and q$ = qo > 2. Thus £^ 2 ^ and S^ 1 '* , vanish in ([4U]) and the asymptotic behavior of 

results from £^ ■ + X„ ■ and X^ 3 -* ■ given in Proposition 16.21 and 16.61 respectively. These 
propositions apply because we assume ([31]) and M > if + <5(go) in Theorem 13.11 Now the 
ratio of the convergence rates appearing in these propositions reads 

l/2-d -(S(q eo )+5(g io +l)+K) ( - u \ 1 / 2 - d 



Hence Case (a) of Theorem 13.11 corresponds to 

and Case |(b)| to 

E (0) ■ + ■ = 0P ( S (3) -1 . 
The proof of Theorem 13.11 follows. □ 

7.2. Proof of Theorems [32] and \3M Here Condition (JSgJ) holds, so that q = 1, q x < 
1/(1 — 2d) and £q = 00 (or equivalently / is an empty set). In particular E„ • vanishes in (|38"|) 



and the asymptotic behavior of S n , j is obtained from those of si 1 '^, ■ and E„ „-. 

Since moreover M > K + d, Proposition 15.21 applies . Using the definition of q$ in (j60p we have 

= qi, and since M>K + d>K + 5{q^) Propositions 16.21 also applies. Finally, observing 
that here too defined in (|26p equals 1 and that M > K + d, Proposition 16.41 applies. Thus, 
using (|46p . it only remains to compare the convergence rates in these propositions. 

We first prove Theorem 13. 2[ Recall that, by Lemma [3TTT since q\ < q*, one has 

v\ < v 2 < u 3 , 

where these three indices are defined in (|38p . In Figure [H we provide pairwise comparisons of 
the rates of convergence of si 1 '^, X„ „• + X„ „• and xi 2 '' •. We obtain domains separated by 
the three curves to = j^ 1 , n~ = 7^ 2 and to = 7^ 3 . Each curve is concerned with a pair of two 
terms among the three and separates the plane (7^, to,-) in two domains, where one of the two 
terms dominates the other. We indicated the dominating term by G for the asymptotically 
Gaussian term S^' 1 -' ^, R for the asymptotically Rosenblatt term „• + E„ A and H for the 
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7 



II, 




Figure 1. Pairwise comparisons of the rates of convergence of S 



,(o) 



^n-,j + S Jy j ( R ) and ^n-,j ( H ) in the P lane H VerSUS n 3 



,(2) 



(1,1,0) 



(G), 



(3) 

term S^ ■ belonging asymptotically to a chaos of order greater than 2. We get three domains 
where one term dominates over the other two: 



7J 3 <C Hj-: since the domain lies both on the left-hand side of the curve rij = 7J 3 and 



on the left-hand side of the curve rij = t^ 2 , R dominates H and R dominates G, hence 
R dominates H and G. That is, the two terms si 1 '* and ■ are both negligible 



with respect to S^^ + S^ 1 ^-. By Proposition I6.2| we obtain Case (a) of Theorem 



,(i) 



jj 1 -C rij <C 7 ? - 3 : since the domain lies both on the right-hand side of the curve 

7J 1 , H dominates R and 

H dominates G, hence R dominates R and G. That is, the two terms S„\ and 



n j 



7^ 3 and on the left-hand side of the curve rij 



+ S„ j are both negligible with respect to S^^-. ^ Proposition 16.41 we obtain 



E (o) 

nj,j 

Case |(b)| of Theorem 



(2) 

3 ' 



n, <C 7^: G dominates -ff and i?, that is, the two terms S^ ** • + Si,' 1 *' „• and Si," 
are both negligible with respect to S, 



,(i) 



,(2) 



(1,1,0) 
nj,3 ' 



By Proposition 15.21 we obtain Case (c) of 



Theorem [ 

This completes the proof of Theorem 13.21 

The proof of Theorem 13.31 is similar except that the assumption qi > q\ implies that 

^3 < ^2 < v\ ■ 

The domains of convergence are now obtained from Figure [2j □ 



8. Technical lemmas 
The following lemma is used in the proof of Proposition 15. II and in that of Lemma 18.31 
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Figure 2. Domains of convergence for Theorem 13.31 



Lemma 8.1. Define, for all a > and Pi G (0, 1), 

■7i,aO»i;0i) = M~ ft , SiGR, (102) 
and, /or any integer p > 2 and /3 = • • • , /3 p ) G (0, l) p , 

Mp—l)a ra P 

Jp,a(si\P)= •••/ [J|s J _i-s l |- ft - 1 |s p |-^ ds p ...ds 2 , si€R. (103) 

J s 2 =— (p— l)a J s p =— a j_2 

T/ien 

(%) «/ /3i + • • • + /3 P > p — 1, we /iaue 



C7 p (/3):=sup sup (\ ai \-(P-^+"'+M) J Pja ( Sl ; /?)) < 

a>0 |*l|<JJO V 7 



00 



f«j if /3i + • • • + f3 p = p — 1, we have 
C p (f3) := sup sup 



-«7p,a(si;/3) < OO , 



a>0 |si|<pa 

Vl + log(pa/|si|) 

(Hi) if there exists q G {2, . . . ,p — 1} snc/i i/iai /3 g + • • • + j3 p = p — q, we have 



C P (P) := sup sup 

a>0 |si|<po V 1 + lQ : 



;(pa/|si|) 



00 



(iv) if Pi + ■ ■ ■ + P p < p — 1 and /or aZZ 5 G {1, . . . ,p — 1}, we have P q + ■ ■ ■ + P p / p — q, 



we have 



C7 p (/3):=sup sup ( a -(P-i-(/Ji+-+0 P )) J p _ a ( Sl ; /?)) < 

a>0 |si|<po V 7 



OO . 
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Moreover, in the case where all the components of f3 are equal to b 6 (0,1), there exists a 
constant c > depending only on b such that 

supc- p (p\) b - 1 C p (b\ p ) < oo , (104) 
p>i 

where l p denotes the p-dimensional vector with all entries equal to 1. 



Remark 8.1. As in lClausel et al.l [2011bJ, all the cases can be compactly written as 

/ -(p-i-G9i+-+A,))+| Si |(p-i-G9 1 +-+A,))- \ 

C p (P) =sup sup — — — - — T — -— J p ,a{si\P) 

a>0| si |< pa \ (l + log(pa/| Sl |)) £ J 

where e = 1 if there exists q G {1, • • • ,p} such that f3 q + • • • + j3 p = p — q and e = otherwise, 
and x+ = max(x,0), x_ = max(— x,0). Now, observing that 

(p - 1 - 2pd) + = (p(l - 2d) - 1)+ = (-28(p))+ = 25„(p) , 

and, similarly, (p — 1 — 2pd)- = 25 + {p), Inequality f j 1 04 [) with b = 2d £ (0, 1) implies there 
exists a constant c > depending only on d such that for any a > 0, |si| < pa 

J p>a ( Sl ; 2dl p ) < C P(p!) 1 - M a 2<5 -W| Sl |- 2<5 +W(l + log(pa/| Sl |)) £ ^ , (105) 

where e(p) is here defined by (|66p . which corresponds to the e above in the case /3\ = • • • = 
/3 P = 2d. 

Proof. Observe first that for all p > 1, 

Jp,a( s r,P) = } l s 2 - sil"^ 1 J p _i j0 (s 2 ;^) ds 2 , (106) 

•/ S2=— (p— l)o 

where /3' = (/32, • • • , /3 P ). The finiteness of the bounds C p (f3) for any integ er p and any 



/3 £ (0, l) p is then proved by induction on p in the different cases in Lemma 9.3 of lClausel et al 
" 2011b( | 



Finally we show the uniform bound (|104[) . that is, that C p (b, . . . ,b) = 0{c\{p\) 1 ~ h ) as 
p — > oo for any fixed b £ (0, 1). We provide a proof only in the case where 1/(1 — b) is not 
an integer (to avoid cases (ii) and (iii)). The proof is similar in the other case. Hence we 
use the induction step described in Case 1 above. Observe that there exists some integer po 
depending only on b, such that for any p > po we have (p — l)b < p — 2, which corresponds 
above to /?2 + • • • + Pp < V ~ 2 (case (iv)). Hence using the induction assumption (|106p . the 
finiteness of C p in case (iv) and the fact that < pa, we get that there exists some positive 
constant c depending only on b such that, 

J p ,a(si;b,--- ,6) < C p _i(6,--- ,b)aP- 2 -^ b ( [ iP \s 2 - Sl \- b ds 2 ) 



< Cp_i(6, •••,&) aP- 2 -^- 1 ^ x (c(b){(2p - l)a)^ b 
= (cCby^Cp-iCb,...^)) 



a P-i~ P b 



This yields that for any p > po(b), C p {b, . . . , b) < c(b)p l b C p -i(b, . . . ,b). Since this holds for 
any p > po(b), the bound f j 1 04 [) follows by induction. □ 



The following lemma provides bounds of a and (3 defined in (|64|) and (J65J). It is used in 
the proofs of Propositions 16.11 16.31 and 16.51 
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Lemma 8.2. One has 



(1) Assume that d > 1/4. Then for any (g, g') E 



•j2 



inf (a{q,q',p)) > 1 - 2d , (107) 

0<p<min(gVg'-2,gAg') 

In any case, 

a(g,g',min(gVg'- l,g Ag')) > 1/2 - d . (108) 

(%) For any g E N 

inf (a(g,g,p)) > min(2(l - 2d), 1/2) . (109) 

0<p<g— 2 

Further, 

a(q,q,q- 1) = min(l -2d, 1/2) . (110) 

f5) For any g E N 



inf (a(g + l,g,p)) > min(3/2(l - 2d), 1/2) . (Ill) 

0<p<q— 1 



C^J For any g E N 



sup (f3(q,p)) < 6+(q) . (112) 
0<P<<? 



Proof. (1) Let us fix (g,g') E N 2 and assume that g' < g. Since the map 

m i— )■ 5+ (m) = max(dm — (m — l)/2, 0) , 

is non-increasing with range in [0, 1/2], one has for < p < min(g — 2, g') 

a(q, q',p) = min(l - 5+(q - p) - 5+(q' - p), 1/2) > min(l - 5 + (2) - 1/2, 1/2) . 

If d > 1/4, 5 + (2) = 2d - 1/2 and thus 

a(q,q',p) > min(l - 2d, 1/2) = 1 - 2d , 

which proves ()107p . Finally, if p = q — 1 and p < q', 

a{q, q',p) = min(l - 5 + (q - p) - 5+(q' - p), 1/2) > min(l - 5 + {l) - 1/2, 1/2) 

= min(l/2 - d, 1/2) = 1/2 - d , 

which proves f|108[) . 
(2) Let us fix g E N, then for any p < q — 2, 

a(q,q,p)=mm(l-6 + (q-p)-S + (q-p),l/2)>min(l-2S + (2),l/2) . 

If d < 1/4, 5 + (2) = and we get a(q,q,p) > 1/2 > min(2(l - 2d), 1/2). If d > 1/4, 
25+ (2) = 25(2) = 4d- 1 and 

a(q,q,p) > min(l - (4d - 1), 1/2) = min(2(l - 2d), 1/2) , 

which gives ()109|) . To prove (|110p , we observe that if p = q — 1, 

a(q,q,p) =min(l-2<5 + (l),l/2) = min(l - 2d, 1/2) . 
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(3) Let us fix q 6 N, then for any p < q — 1, 

a(q + 1, q,p) = min(l - 6+(q + 1 -p) - 5+(q - p), 1/2) > min(l - <5+(2) - S+(l), 1/2) 

= min(l-d-5+(2),l/2) . 

If d < 1/4, <5 + (2) = and a(g + l,q,p) > min(l - (2,1/2) = 1/2. If d > 1/4, 
<5+ (2) = 2d - 1/2 and (jTTTI) follows from 

a(g + l,g,p) > min(l - d - (2d - 1/2), 1/2) = min(3(l - 2d)/2, 1/2) . 

(4) If [3(q,p) = 0, then (3(q,p) < 8+(q). Now consider the case where 

/% p) = max(5 + (p) + 6+(q - p) - 1/2, 0) > , 

that is, 5+(p) + 5+(q — p) — 1/2 > 0. In this case, 5+(p) and 5+(q — p) are both positive 
(since < <5+(-) < 1/2) and they respectively equal 5(p) and S(q — p). Then we obtain 

max(Mp) + M<7 " P) ~ 1/2, 0) = S(p) + % - p) - 1/2 = 5(q) , 
which again implies (jll2D , 

□ 

The following result provides a bound of def ined in (|53[) , in the case where p > 0. It is a 
refinement of Lemma 10.1 of Clausel et al. 2011bl |. It is used in the proof of Proposition 15.31 

Lemma 8.3. Suppose that Assumptions A hold and let p be a positive integer. Then there 
exists some C > neither depending on p nor j such that for any (6,6) £ ^ 2 > 
(i) if for any s £ {1, • • • ,p}, s(l — 2d) ^ 1 then, 

2(S+(p)+K) 

| S rte. & )| S ^)- (1+7;|Ki}|) l M(1+7;|te}|)Mp) . am 

(ii) if there exists s £ {1, • • • ,p} such that s(l — 2d) = 1, then, 

|«f (6,6)1 < ^(p!) 1 - 2 " 7** MTi) • (H4) 



Remark 8.2. In Case (ii) of Lemma 18.31 we have p > 1/(1 — 2d), hence <5+(p) = 0. Equa- 
tions (|114p and (|1 13j) can thus be written as a single bound, namely, 

2(S+(p)+K) 

tf>te,6)i < ^ (r^iFWT^TfeTF^ (log ^ w ' (115) 

where e(p) is defined by ([66j) . 

Proof. By (27r)-periodicity of re^ (6,6) along both variables 6 and 6, we ma y take 6,6 £ 
[— vr,7r]. The remainder of the proof shows that f j 1 1 5 j) holds for such (6,6)- 
Note that by assumption, 

/(A)<C|A|" M , 
where C > only depends on /*. Using (fl37D . (I53]l and (fT03l) with 

Mi = H h Ap) , 

we get 

P U < r^ 2 ^+ 5 ^) ^r(Pjj2dlp)d^ 

•/-pis* Il<=i (1 + Ii Kmi/Tj + 611) 
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Then, by (|105p . there exists C > not depending on j,p such that, for all (^1,^2) G [ — 7r 7 7r ] 2 i 

-AM < r*r n ni-W 2K+2( 5 ( P ) +fi _(p)) |/ii|- 2 ^ (p) (l + log(P7WI/*i|)) e(p) dm 

J-pj^ Ui=i (1 + 7j I W7j + 
Using that 6{p) = 5+(p) — 5-{p) and the Cauchy-Schwartz inequality, to obtain (|115p . it is 
sufficient to show that, for all £ G (— it, 7r], 

where C is a positive constant. 

If S(p) > the rest of the proof is similar to that of Lemma 10.1 in Clausel et al. 2011bl ] 
and is thus omitted. 

We now take 5+(p) = 0, so that (|116p becomes 

The denominator in the integral is a (27T7j)-periodic function of fix, hence the integral over 
[~Plj' K iPlj'K\ is bounded by the sum of at most p + 1 integral of the form 

l + log^-TT/M^dMl 



+ 2/c7j7r) with I(y) = j 

J M 



, A(y) (l + | W - y |)2(^+ Q ) ' 

where k £ Z and A(y) = [— pjjir,p~fjir] n (y — 7j7r, y + 7j7t]. We observe that I(y) is maximal 
at y = where it takes value 



Since the last integral in the previous display is finite for e(p) = 0, 1, we finally obtain (|117p . 

□ 

Appendix A. Integral representations 

It is convenient to use an integr al representat i on in the spect ral domain to represent the 
random processes (see for example Maiorl 1981 ]. Nualart 20061 ]). The stationary Gaussian 
process {X^, k G Z} with spectral density ([2]) can be written as 

e iX£ f 1 / 2 (X)dW(X)= f - iX J <iW(X), £gN. (118) 

-7T J— TT \± e I 

This is a special case of 

T(g) = [ g(x)dW(x), (119) 

where W(-) is a complex- valued Gaussian random measure satisfying, for any Borel sets A 
and B in R, E(W(A)) = 0, E(W{A)W(B)) = \AnB\ and 



W(A) = W(-A) . 

The integral (|119p is defined for any function y G L 2 (M) and one has the isometry 

n\i(g)\ 2 ) = I \g(x)\ 2 dx . 
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The integral 1(g), moreover, is real-valued if 



a(x) = g(-x) . 

We shall also consider multiple ltd- Wiener integrals 

T q (g)= [ g(X 1 ,--- ,X q )dW(X 1 )---dW(X q ) 

where the double prime indicates that one does not integrate on hyper diagonals Aj = ±Xj,i ^ 
j. The integrals I q (g) are handy because we will be able to expand our non-linear functions 
G(Xk) introduced in Section Q] in multiple integrals of this type. 

These multiples integrals are defined for g G L 2 (K 9 ,C), the space of complex valued func- 
tions defined on M q satisfying 



g(-xi, ■■■ , -x q ) = g(x 1 , ■ ■ ■ , x q ) for (x 1} ■ ■ ■ , x q ) G W , (120) 
IMli 2 := / \d(xi,--- ,x q )\ 2 dxi- ■ -dx q < oo . (121) 

Hermite polynomials are related to multiple integrals as follows : if X = J K g(x)dW(x) with 
E(X 2 ) = J R g 2 (x)dx = 1 and g(x) = g(—x) so that X has unit variance and is real-valued, 
then 



H q (X) = I q (g* 



g(x 1 )---g(x q )dW(x 1 )---dW(x q ) 



(122) 



Since X has unit variance, one has for any i G 



H q (X e ) — H q 



J& f 1/2 



(C)dW(0 



'(-• 7T,7r]8 

Then by ([40 p . we have 



x ( / 1/2 (£i) x ••• x ^(Q ) dWfo) ■ ■ ■ dW(Q . 



(9)> 
k' 



(123) 



with 



M 
because 



• • • = e i ^« 1+ - + ^) x hf >(& + • • • + e 9 )/ 1/2 (6) • • • Z 1/2 fe)l^)(e) , (124) 



Y^^ + -^Kf\ l3 k-i) 



(el 



^7jfc(ei+-+€,)^ e -i«(€i+"-K,) h W( 

e i 7 ,fe(a+-+c 3 )hW(^ 1 + ... + ^ 



by(H2ZD. 

The following proposition can be found in lPeccati and Taaaul [20111 ]. F ormula (9.7.32) . It 
is an extension to our complex- valued setting of a corresponding result in iNualart 2006 ] for 
multiple integrals in a real-valued setting. 
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Proposition A.l. Let (q,q') G N 2 . Assume that f,g are two symmetric functions belonging 
respectively to L 2 (R q ) and L 2 (R q ') then the following product formula holds : 

q 



Uf)iAa) = (JJ (JJ i q+q >- 2p (f® P g), (125) 

where for any p G {1, • • • , g A q'} 
{f® P 9)(h, ■■■ , t q+q >-2p) = (2vr) p / /(ii, • • • , s)5(tg_ p+ i, • • • , t q+q ,_ 2p , -s)d p s . (126) 

Appendix B. The wavelet filters 

The sequence {Y t } t£ % can be formally expressed as 

Y t = A~ K G(X t ), t€Z. 

The study of the asymptotic behavior of the scalogram of {Y t }tez at different scales involve 
multidimensional wavelets coefficients of {G(X t )}tez and of {Y t } t ez- To obtain them, one 
applies a multidimensional linear filter hj(r),r G Z = (hj^(r)), at each scale j > 0. We shall 
characterize below the multidimensional filters hj(r) by their discrete Fourier transform : 

n,(A) = ^h J (r)e- iA -, A G [-tt.tt] , h,(r) = -L f h,(A)e iA MA,r G Z . (127) 

The resulting wavelet coefficients W^fc, where j is the scale and k the location are defined as 

W jjk = £ h ihjk - t)Y t = Y, hji-yjk - t)A- K G(X t ), j>0,k£Z, (128) 
tez tei 

where 7 j t oo as j t oo is a sequence of non-negative scale factors applied at scale j, for 
example 7 j = 2 J . We do not assume that the wavelet coefficients are orthogonal nor that 
they are generated by a multiresolution analysis. Our assumption on the filters h, = (hj^) 
are as follows : 

(W-a) Finite support: For each £ and j, {hj t e(T)} T £z has finite support. 
(W-b) Uniform smoothness : There exists M > K, a > 1/2 and C > such that for all j > 
and A G [— it, tt], 

C 7l 1/2 | 7j A| M 

|h ' (A)l £ (i ^iAlP= ' (129) 
By 27r-periodicity of /ij this inequality can be extended to A G E as 



7 J - /2 |7 J -{A}| M 
(1 + 7 ,|{A}|)^ ' 



h *wi < g ^ . , n +M ■ ( i3 °) 



where {A} denotes the element of (—tt, tt] such that A — {A} G 2ttZ. 
(W-c) Asymptotic behavior: There exists a sequence of phase functions $j : E — > (—tt, tt] 

and some non identically zero function hoo such that 

.lim ( 7 7 1 /%( 7 7iA)) = n 00 (A), (131) 

j — >+oo J J 

locally uniformly on A G E. 
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In (W-c) locally uniformly means that for all compact K C 



sup 



-1/2 



Assumptions (|129j) and (|13ip imply that for any Agl, 

|A| M 



|hoo(A)| < C- 



l + |A|) a + M 



(132) 



Hence has entries in L 2 (R) . We let be the vector of L 2 (M) inverse Fourier transforms 
of hi oo , that is 



= d(hoo)(0 = [ h 00 (t)e' i ' T « d% £ G 



(133) 



is defined for any f £ L 2 

Observe that while hj is 27r-periodic, the function hoc has non-periodic entries on 
For the connection between these assumptions on hj and corresponding assu mptions on the 



scaling function ip and the mother wavelet if) in the classical wavelet setting see lMoulines et al 



2007I ] . In particular, in that case, one has h^ = ip(0)ip. 



A more convenient way to express the wavelet coefficients Wj & defined in (|128p is to 
incorporate the linear filter A - ^ in (|128p into the filter hj and denote the resulting filter 



h[ K) . Then 



where 



hf\\) = (1 - e- iX )-%W 
IK) 



is the discrete Fourier transform of h^- . Using (|130p we get 

u 2+K | 7 i{A}| M -* 



£f°(A) 



<C 7 



A G R, j > 1 



(1 + 7 ,|{A}|)^' 
In particular, since we assume if M > K, we get 

hf } (A) <C^ /2+K (l + ^\{X}\)- a - K , AGM,j>l 



(134) 



(135) 



(136) 



(137) 



By Assumption (|129p . hj has null moments up to order M — 1, that is, for any m G 
{(),•••, M-l}, 



(138) 



Observe that A^y is centered by definition. However, by (I138p . the definition of W^fc only 
depends on A M Y. In particular, provided that M > K + 1, its value is not modified if a 
constant is added to A K Y, whenever M > K + 1. 
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